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The anisotropic optical properties of multiferroic BiFeOs thin films have been determined with 
Mueller matrix ellipsometry at room temperature. The full dielectric tensors of tetragonal-like and 
rhombohedral-like BiFeOs phases epitaxially grown on LaAlOs and SrTi 03 single crystal substrates, 
respectively, within the spectral range of 0.6 and 6.5 eV are reported. Strain-driven anisotropy 
changes and transition shifts are observed as well as evidence of sub-band gap many-particle excita¬ 
tions are found. The transition shifts, mostly to higher energies for the highly strained tetragonal-like 
BiFeOs phase on LaAlOs, are indicative of band structure differences. Additionally, optical mod¬ 
eling, confirmed by piezoelectric force microscopy studies, revealed that the average polarization 
direction of bivariant BiFeOs on LaAlOs is not parallel to the crystallographic [001] direction but 
tilted by about 7°. Spectral weight analyses reveal phase-dependent differences, underlining that 
theoretical calculations of optical spectra need further improvement to appropriately account for 
electronic and excitonic correlations to fully understand multiferroic BiFeOs- 

PACS numbers: 78.20.Ci, 78.20.Fm, 07.60.Fs, 75.85.+t, 81.05.Xj 


I. INTRODUCTION 

Single crystalline bismuth ferrite (BiFeOs) is a mul¬ 
tiferroic perovskite structure and exhibits magnetic as 
well as strong ferroelectric behavior at room tempera¬ 
ture. For about the past decade, BiFe0 3 has been of 
strong research interest due to its potential applicability 
in ferroelectric memory devices and spintronics as well 
as photovoltaics, for example^. While the lattice system 
of bulk BiFe0 3 is rhombohedral, the crystal structure 
of thin films can be engineered by introducing epitaxial 
strain. Depending on the choice of single crystalline sub¬ 
strate materials and their different lattice parameters, the 
thin film BiFe0 3 crystal structure and associated physi¬ 
cal properties, such as transition energies, can be modi¬ 
fied^—. 

Although there are already numerous publications 
about the optical properties of bulk or thin film BiFe0 3 , 
often the highly anisotropic nature of crystalline BiFe0 3 
has not been appropriately considered during experiment 
and data analysis. Apart from that, especially the optical 
band gap energy and its absorption onset, and whether 
BiFe0 3 is a direct or indirect material is controversially 
discussed. 

Most ab initio calculations agree that BiFe0 3 is a semi¬ 
conductor with an indirect band gap that is very close 
to the first direct transition due to the flatness of the 
bands^. However, while some argue that the closeness of 
direct and indirect gaps is due to the valence band being 
very flat^i, others have calculated very flat conduction 
bands^. In general, ab initio calculations to character¬ 
ize multiferroics with transition metal cations and oxy¬ 
gen may be very complicated as both the exchange in¬ 
teraction and electron correlations have to be taken into 
account. Additionally, excitons are typically not consid¬ 
ered when calculating optical properties, even though ex¬ 


citonic effects may strongly influence the dielectric func¬ 
tion tensor—. Given the multitude of available variations 
of computational modeling with density functional theory 
(DFT), different results are not surprising^. 

Two experimental reports regarding polycrystalline 
BiFe0 3 thin films suggest the presence of an indirect gap 
roughly 1.0 eV below the first direct transitio n 12 ! 13 . How¬ 
ever, many others argued that no indications of an indi¬ 
rect gap were found, and concluded that BiFe0 3 is a di¬ 
rect band gap material with a transition energy between 
2.6 and 3.1 eV at room temperature (mainly depending 
on the crystalline phase)^ - — . 

Interestingly, all of these aforementioned references 
point out that the absorption onset starts significantly 
below the first direct transition. One of the reasons for 
this can be a shallow oxygen vacancy state below the con¬ 
duction band, as calculated by Clark and Robertson^ 
and as is consistent with observations by Hauser et ali^. 
Ju and Cai showed in a theoretical study that the absorp¬ 
tion onset might strongly redshift with increasing defect 
states in the form of oxygen vacancies^!. However, even 
though a redshift of the first allowed direct transition was 
confirmed with increasing oxygen vacancies a shift of the 
absorption onset could not be observed experimentally^. 
The early absorption onset far below the main transition 
(often observed as a long structureless enhanced spectral 
weight) was assigned by Pisarev et al. to charge transfer 
instabilities accompanied by a self-trapping of excitons^. 
Xu et al ., however, observed two distinct peaks at around 
1.4 and 1.9 eV in transmittance measurements of a bulk- 
like single crystal and assigned them to on-site crystal- 
field transitions^. These peaks are consistent with pre¬ 
viously reported many-particle transition bands also in¬ 
volving excitons^. 

It is noteworthy that in most previous studies the op¬ 
tical properties have been determined by making use 
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FIG. 1. In-plane PFM images of BiFe03 on LaA103 (left) 
and SrTiOs (right). The schematic insets in the top left cor¬ 
ner depict the polarization directions. The insets in the top 
right corner show the respective single-domain out-of-plane 
images with indicated polarization directions. The scale bars 
are 500 nm. 


of parametrized oscillator models to extract the dielec¬ 
tric function from spectroscopic ellipsometry or trans¬ 
mittance data, for example^ " 22 ! 25 ! 26 . Such sub-band gap 
transitions have not been reported in these studies. In 
general, with the use of parametric physical line-shape 
models, a certain risk is involved for subtle spectral fea¬ 
tures to be neglected by the line shape of the model func¬ 
tion. Since the characteristics of indirect transitions are 
often slight and the absorption due to crystal-field transi¬ 
tions is approximately three orders of magnitude smaller 
than the absorption above the charge gap, particular care 
must be taken here^i. 

Here, we present the anisotropic optical properties of 
multiferroic BiFeOs thin films with a nominal thickness 
of around 35 nm as determined with Mueller matrix ellip¬ 
sometry at room temperature. Ellipsometry within the 
Mueller-Stokes formalism has been shown to be an excel¬ 
lent technique for the determination of the dielectric func¬ 
tion tensor of biaxially anisotropic materials^ - —. The 
full dielectric tensors as well as major polarizability direc¬ 
tions of epitaxial rhombohedral-like and tetragonal-like 
BiFeOs phases are discussed. Charge transfer transitions 
of both films are quantified and strain-induced differences 
in excitation energies between rhombohedral and highly 
strained tetragonal phases are reported. A comparison of 
experimental data from both crystalline phases with re¬ 
cent DFT with the Heyd-Scuseria-Ernzerhof hybrid func¬ 
tional is made to reveal that many-body effects are not 
yet appropriately accounted for and particularly excita¬ 
tions of an excitonic nature need to be taken into account 
when calculating the dielectric function tensor—. 


II. SAMPLE DETAILS 

Epitaxial BiFeOs thin films were deposited on LaAlOs 
(001) and SrTiOs (001) single crystal substrates by 
pulsed laser deposition. While the LaAlOs substrate has 
no intentional miscut, SrTiOs has an intentional 4° mis- 
cut along the [110] c direction. A stoichiometric target 
was ablated by using a KrF excimer laser (A = 248 nm) 
with an energy density of ~ 1.5 J/cm 2 and a repetition 
rate of 10 Hz. The deposition temperature was 650 °C 


and the oxygen partial pressure was 100 mTorr—. 

Figure [T| depicts in-plane and out-of-plane (insets) 
piezoelectric force microscopy (PFM) images of BiFeOs 
films on LaAlOs and SrTiOs, respectively. The left image 
(BiFeOs on LaAlOs) shows part of one bivariant in-plane 
domain and the diagonal “stripes” are parallel to atom¬ 
ically flat LaAlOs terraces. Due to the lattice constant 
of LaAlOs (a = 3.79 A), the coherent BiFeOs him is 
in a highly strained tetragonal-like monoclinic phase (T 
phase) with a giant c/a ratio (1.23) 37 39 . Based on the 
terrace widths of approximately 76 nm an unintentional 
miscut angle of 0.28° along the [110] pc direction can be 
calculated. The BiFeOs him on the SrTiOs substrate is 
not as strained and is in a rhombohedral-like monoclinic 
phase (R phase) with a c/a ratio of 1.03^. The thin film 
is partially relaxed due to the large substrate miscut an¬ 
gle. A near perfect single-domain character with approx¬ 
imately 5% other random domain variants is observed in 
large scale in-plane PFM images. For both rhombohedral 
and tetragonal BiFeOs samples, a single-domain out-of¬ 
plane character is confirmed. 

The monoclinic distortion in both cases is very small 
(< 2°) and therefore an orthorhombic lattice has been 
assumed for the following optical analysis^. 


III. SPECTROSCOPIC ELLIPSOMETRY 

Spectroscopic Mueller matrix ellipsometry spectra 
within the spectral range from 0.6 to 6.5 eV in steps of 
20 meV were acquired using a commercial rotating ana¬ 
lyzer instrument with a compensator (VASE, J. A. Wool- 
lam). The samples were mounted on a precision rotation 
stage (RS40, Newport) to perform azimuth-dependent 
measurements, and the in-plane rotation angle </> was var¬ 
ied from 0° to 320° in steps of 40°. At each in-plane 
orientation, data were taken at three angles of incidence 
4> a (50°, 60°, 70°). Such an angle-resolved measurement 
scheme is necessary for a complete characterization of ar¬ 
bitrary optically anisotropic samples^ - — . Ellipsometric 
spectra (4/ and A) for the pristine isotropic substrates 
were measured in the same energy range at a single in¬ 
plane orientation ^ 41 i 42 . 

The complex dielectric function for the isotropic 
SrTiOs substrate has been calculated by wavelength- 
by-wavelength inversion of the experimental data (pseu¬ 
dodielectric function)^. For LaAlOs the pseudodielec¬ 
tric function has been further parameterized by a physi¬ 
cal line-shape model in order to avoid experimental data 
noise from becoming part of the calculated dielectric func¬ 
tion and thus any further analysis of the BiFeOs thin 
fil m 41 ! 43 . 

The stratified layer optical model for both samples un¬ 
der investigation comprises, besides the substrate, two 
additional layers accounting for the BiFeOs thin film and 
the surface roughness, respectively. The optical model 
equivalent of surface roughness is a layer with thickness 
d and optical constants typically calculated based on a 
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FIG. 2. Experimental (dotted lines) and best-match calcu¬ 
lated (solid lines) normalized Mueller matrix spectra Mij /Mu 
for BiFeOs on LaA 103 (top) and BiFeOs on SrTi 03 (bot¬ 
tom). The data are plotted for an angle of incidence = 70° 
and in-plane orientations of p — 20° and 40° for the LaAlC >3 
and SrTiOs samples, respectively. Note that the off-diagonal 
Mueller matrix elements M13, M14, M23, M24, M31, and M32 
are plotted xlO for clarity. 


Bruggeman effective medium approximation with 50% 
material and 50% void (Rej^} = 1 and Imj^} = 0)42. 
The BiFeOs thin film in both cases is assumed to have 
an orthorhombic crystal system with complex dielectric 
function parameters e a , £&, and e c along the major polar¬ 
izability axes a, 6, and c, respectively. Real-valued and 
wavelength-independent Euler angles (/?, 0, and ^ are used 
to rotate between the Cartesian laboratory frame and 
the Cartesian polarizability system and, together with 
a layer thickness d, complete the list of model parame¬ 
ters^ - —. Here, the complex dielectric function parame¬ 
ters have been calculated wavelength by wavelength, i.e., 
£ a , and £ c are obtained at each wavelength, indepen¬ 
dent from all other spectral data points. Hence, no phys¬ 
ical line-shape model was implemented, which allows for 
an unbiased extraction of the intrinsic dielectric function 
tensor. Independent Kramers-Kronig consistency tests 
can then be done with the individual dielectric function 
parameters^. 

Data analysis requires nonlinear regression methods, 
where measured and calculated Mueller matrices are 



FIG. 3. Real and imaginary parts of the dielectric func¬ 
tion tensor obtained from best-match model calculations for 
BiFeOs on LaAlOs (top) and BiFeOs on SrTiOs (bottom). 
The insets show schematically the orientation of the major 
polarizability axes a (black), b (red), and c (green). 


matched as closely as possible by varying the above- 
mentioned model parameters, thereby minimizing a 
weighted test function. Both samples discussed here have 
been analyzed using a multisample-configuration analy¬ 
sis scheme, i.e., all nine in-plane orientations have been 
included in the regression analysis and only the azimuth 
Euler angle ip has been set according to the measurement 
configuration {<p n +i = p n ± 40°)2i. 


IV. RESULTS AND DISCUSSIONS 

Figure [2] depicts representative experimental and best- 
match calculated Mueller matrix spectra for both samples 
under investigation. Spectra are plotted for one angle 
of incidence and one in-plane orientation (see the figure 
caption for details). Both samples exhibit anisotropic op¬ 
tical behavior, as is evident from the off-diagonal Mueller 
matrix elements, and the degree of anisotropy is larger 
for the rhombohedral BiFeOs film. As a result of the 
nonlinear regression fitting, film thicknesses of 40.8 ±0.1 
and 32.7 ± 0.3 nm have been determined for BiFeOs on 
SrTiOs and LaAlOs, respectively. The latter is in very 
good agreement with a thickness of 31 nm measured by x- 






























4 


ray diffraction. The thicknesses of the surface roughness 
layers are 2.5 db 0.1 and 3.0 db 0.2 nm, respectively. 

A. Dielectric function tensors 

The wavelength-by-wavelength extracted dielectric 
function tensor as well as a schematic representation of 
the major polarizability coordinate system orientation 
with respect to crystallographic axes are presented in 
Fig. [3l In both cases a large degree of birefringence and 
dichroism is observed, and while BiFeOs on LaAlOs in¬ 
dicates a uniaxial character (e a « % s c ), BiFeOs on 

SrTiOs has biaxial optical properties (e a / £5 / e c ). 
Within the investigated spectral range BiFeOs cm SrTiOs 
exhibits an average of 30% and 44% less birefringence and 
dichroism, respectively, compared to the highly strained 
BiFeOs on LaAlOs. Similar to previous reports on bulk¬ 
like rhombohedral BiFeO s 20 ! 45 , a negative birefringence 
(n c < n a ,nb, rij = Rely^}) is observed below 2.5 eV 
for both films. Interestingly, the birefringence is sig¬ 
nificantly more pronounced for tetragonal-like BiFeOs 
(~ 70% at 1.5 eV), which is in contrast to first-principles 
DFT, where the birefringence was calculated to be nearly 
identical between both phases^ While the overall shapes 
of both dielectric function tensors show similarities, the 
peak positions in the case of BiFeOs on LaAlOs are 
blueshifted with respect to BiFeOs on SrTiOs. Besides 
that, the behavior in the low-energy range shows some in¬ 
teresting characteristics: On LaAlOs, Im{£ a } and Imje^} 
below 1.8 eV and Im{s c } below 1.2 eV are zero (within 
the experimental error), while this is not observed within 
the measured spectral range for the film on SrTiOs. 

Furthermore, optical modeling revealed that even 
though BiFeOs on LaAlOs exhibits uniaxial-like optical 
properties, the major polarizability axis e c (optic axis) 
does not coincide with the [ 001 ] crystallographic axis 
but rather is tilted away from the surface normal by 
6 = 7.4 =b 0.9°. This tilt has been confirmed with careful 
PFM studies as well as measurements of in- and out-of¬ 
plane polarization component s 46 ! 47 , resulting in a value 
of < 9°. The small unintentional substrate miscut results 
in an overall preferred domain alignment responsible for 
this net polarization tilt. The Euler angle ip = 0° was 
not included in the analysis, resulting in the major po¬ 
larizability axis a being in planed. 

For BiFeOs on SrTiOs, the Euler angles 0 = 42.8° and 
'ip = 0 ° have not been included in the regression analysis, 
and the calculated angle ip = 41.865 =b 0.007° indicates 
that the polarization is along the [ 111 ] directio n 48 ! 49 . 

B. Transitions and Excitations 

The dielectric function tensor enables computation of 
many relevant sample properties as well as the absorp¬ 
tion coefficient a, which in turn allows for the determi¬ 
nation of transition characteristics. Plots of the square 
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FIG. 4. Plots of (ahv) 2 and (ahv) 1 / 2 vs photon energy visu¬ 
alizing direct and indirect transitions, respectively, for BiFeOs 
on LaAlOs (top) and BiFeOs on SrTiOs (bottom). The in¬ 
set shows the absorption coefficient a vs photon energy high¬ 
lighting two subtle excitations for R-phase BiFeOs- No such 
excitations are observed for T-phase BiFeOs- 


and square root of ahv allow for a graphical determi¬ 
nation of direct and indirect transitions, respectively, by 
extrapolating the linear regimes to ahv — 05fi. For an in¬ 
direct band gap at room temperature, two linear branches 
should be observable, corresponding to phonon emission 
and absorption^. 

Figure |4] depicts plots of {ahv) 2 and {ahv) 1 ! 2 versus 
photon energy hv. For BiFeOs on LaAlOs (top panel), 
the first direct transition is observed at 3.07 eV by ex¬ 
trapolating both axes a and 6 ; this energy position is 
in agreement with previous reports for a direct gap 19 ! 22 . 
Since this transition is only to Fe 3 d xy orbitals, axis c is 
not considered for the graphical analysis^. For BiFeOs 
on SrTiOs, the first direct transition is found by extrapo¬ 
lation of all three axes at 2.68 eV (Fig. 0J bottom panel). 
This is in agreement with previous reports for compara¬ 
ble R-phase samples^ - — and corresponds to a redshift 
of around 390 meV with respect to BiFeOs on LaAlOs. 

In addition to the first direct transition, there are sev¬ 
eral linear regimes in the ( ahv ) 1 / 2 plot that can be ex¬ 
trapolated to determine other excitations (Fig. 0]). Par¬ 
ticularly, two linear branches are observed, potentially 
representing phonon absorption (E g — E p = 1.86 eV) and 
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emission (E g + E p = 2.21 eV) in T-phase BiFeOs. This 
would set the indirect band gap E g at 2.035 eV with an 
involved phonon of 175 meV. An analogous analysis for 
BiFeOs on SrTiOs would result in an indirect band gap 
at E g = 1.64 eV with a phonon energy of E p = 410 meV. 
However, the required optical phonon energies, partic¬ 
ularly in the case of T-phase BiFeOs, are too high for 
such a scenario. Nevertheless, since the branch related 
to phonon absorption is usually very subtle it would be 
reasonable to assume that extrapolation of a single linear 
regime results in an indirect gap with involved phonon 
emission E g + E p . 52 This, however, leaves the graphical 
determination of a possible indirect band gap inconclu¬ 
sive. Any of the three determined values would be a can¬ 
didate (with sub-band gap excitations where applicable), 
but none of the theoretical calculations have placed an 
indirect gap energetically so far below the first direct gap 
to make an assumption^. 

The graphically determined transitions at 2.35 and 
2.49 eV for the T and T phase, respectively, are responsi¬ 
ble for small shoulders in the dielectric function and have 
been observed in many previous studies. It was proposed 
that these excitations are likely defect related and, due 
to the presence of a moderately shallow oxygen vacancy 
state, approximately 0.3-0.6 eV below the direct band 
gap 15 i 23 . The next lower-energy excitations (2.05 and 
2.21 eV for the R and T phase, respectively) may then 
be assigned to a dipole-forbidden on-center tig(n) —)> t^g 
charge transfer transition, which has been predicted to be 
redshifted by about 0.8 eV with respect to the respective 
direct gap2&. 

Interestingly, for T-phase BiFeOs, two additional sub¬ 
tle excitations below 2 eV are identified and emphasized 
in the inset of Fig. 01 The peak positions have been de¬ 
termined by a dedicated nonlinear regression analysis of 
the wavelength-by-wavelength extracted ImjT} and are 
at 1.37 and 1.86 eV. These energy positions are in very 
good agreement with previously observed sub-band gap 
transitions as a result of transmittance measurements on 
a bulk-like rhombohedral single crystal. The two very 
weak peaks can be attributed to many-particle transition 
bands and specifically comprising a pure exciton tran¬ 
sition 6 Ai g —)> A T\ g and another 6 Ai g —)> ^T^g on-site 
Fe 3+ crystal-field transition barely allowed by spin-orbit 
couplin g 27 ! 28 ! 53 . Together, these sub-band gap excita¬ 
tions are responsible for the absorption onset for T-phase 
BiFeOs being outside the measured spectral range. 

These transitions are not observed for the highly com¬ 
pressed in-plane axes a and b of T-phase BiFeOs, but 
rather the onset of absorption is at around 1.86 eV. For 
the strained out-of-plane axis, however, the absorption 
onset is significantly redshifted. Unfortunately, due to a 
lower experimental sensitivity to axis c and hence some 
data noise, the subtle peak at around 2.1 eV cannot be 
characterized with necessary certainty, but it could be 
another many-particle excitation^. A probable explana¬ 
tion for the difference in the absorption onset is that the 
charge transfer exciton self-trapping, which is partially 



FIG. 5. Imaginary parts of the dielectric function tensor e c 
for R- and T-phase BiFeOs (left) and normalized effective 
carrier number iV e ff (right) of this work and compared to lit¬ 
erature values from Dong et alM. The two diagram insets 
show schematically the dominant excitations. 

responsible for the shallow absorption tail, is governed 
by lattice strain and is mostly suppressed here2£. 


C. Comparison with theoretical calculations 

The comparison of experimentally obtained imaginary 
parts of the dielectric function along axis c and theoreti¬ 
cally calculated values taken from Dong et alS^ is shown 
in Fig. [5] The dielectric function of each phase below 6 eV 
is dominated by two intense bands (around 3 and 4.5 eV 
for R phase, and 3.5 and 5 eV for T phase), which can be 
assigned to O 2 p z — )► Fe 3 d XZjyz and O 2 p z —>> Fe 3 d z 2 elec¬ 
tron excitation s 26 ! 35 . In general, it is important to note 
that a peak in ImjT} does not necessarily correspond to 
a single interband transition only as multiple direct or 
indirect transitions may be found in the band structure 
with an energy corresponding to or very close to the same 
peak 4 

Experimentally, the two main T-phase peak positions 
are blueshifted with respect to the R phase by 0.3 and 
0.8 eV for low- and high-energy excitations, respectively. 
This blueshift can be explained by local symmetry break¬ 
ing of the FeC>6 octahedral due to the high compressive 
strain. 19,35 The thin and thick arrows in the inset of Fig. [5] 
depict schematically the experimentally found relative 
transition strengths for T- and T-phase BiFeOs. The 
additional T-phase electronic feature centered at 6.2 eV 
may be due to a strongly hybridized majority channel 
O 2 p z + Fe ?>d z 2 —Bi 2 p state excitatio m 19 ! 35 It is as¬ 
sumed that for T-phase BiFeOs such a transition is re¬ 
sponsible for the 0.8 eV redshifted small peak at 5.4 eV. 

For T-phase BiFeOs, the theoretically determined low- 
energy p-d transition strength is slightly underestimated 
while the high-energy one is in excellent agreement with 
the experiment; however, the peak position is blueshifted 
by approximately 0.7 eV. Note that all dielectric func¬ 
tions obtained by DFT had been shifted equally to match 
the location of the low-energy band of T-phase BiFeOs. 
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Regarding the energy positions of the main transitions, 
the situation is the opposite for T-phase BiFeOs. While 
the high-energy excitation peak around 5 eV matches 
very well with the experiment, the low-energy p-d tran¬ 
sition is blue-shifted by about 0.5 eV. Both calculated 
excitation strengths are significantly overestimated com¬ 
pared to the experiment. 

In addition, Fig.[5]shows the normalized effective num¬ 
ber of carriers N e ^ defined as the integration of the opti¬ 
cal conductivity a over the measured frequency ranged. 
Here, the averaged optical conductivity along all major 
polarizability axes is taken into account. Results from 
the DFT calculations show that except for a small range 
between 3 and 4 eV, an identical amount of carriers for 
R- and T-phase BiFeOs are participating in optical tran¬ 
sitions. The experimentally observed situation, however, 
is different and the number of effective carriers is signif¬ 
icantly larger for T-phase BiFeOs within the measured 
spectral range. This means a substantial spectral weight 
transfer to higher energies (> 6.5 eV) occurs for T-phase 
BiFeOs. Consequently, the current first-principles DFT 
calculations with exchange-correlation functionals do not 
fully account for all many-body effects that may influ¬ 
ence the electronic structure. Hence, in order to improve 
the theoretical calculations of optical spectra and thus al¬ 
low for a more comprehensive interpretation of this com¬ 
plex multiferroic compound, many-body effects such as 
exchange and correlation interactions must be further re¬ 
fined and especially excitonic excitations taken into ac¬ 
count. 


V. CONCLUSIONS 

In conclusion, we have carefully determined and 
investigated the anisotropic dielectric function tensor 
of tetragonal-like and rhombohedral-like multiferroic 
BiFeOs phases. Optical modeling of spectroscopic 
Mueller matrix ellipsometry data enabled the quantifica¬ 
tion of strain-driven birefringence and dichroism as well 
as band structure changes in terms of transition energy 
shifts. 

It is found that the birefringence in the low-energy 


range (< 2.5 eV) is significantly larger for tetragonal-like 
BiFeOs, which is expected but is against previous calcu¬ 
lations. The first direct transition has been determined 
for each of the two BiFeOs phases, and they are separated 
by approximately 0.4 eV. Further excitations have been 
quantified, and the difficulty of a graphical identification 
of a potential indirect band gap is discussed. The unbi¬ 
ased wavelength-by-wavelength analysis also allowed for 
a characterization of subtle sub-band gap excitations in 
rhombohedral BiFeOs on SrTiOs, which are not present 
in the same manner in highly strained tetragonal BiFeOs 
on LaAlOs. 

Furthermore, optical modeling and polarization com¬ 
ponent measurements revealed that the major polariza¬ 
tion of BiFeOs ° n LaAlOs within a bivariant in-plane con¬ 
figuration is tilted away from the crystallographic [001] 
axis and the surface normal. 

Comparing our experimental data to recent density 
functional theory calculations has exposed shortcomings 
of the used first-principles algorithms. Highly strained 
tetragonal-like BiFeOs on LaAlOs exhibits a substantial 
spectral weight transfer to higher energies beyond 6.5 eV 
compared to its rhombohedral-like counterpart, which is 
not predicted by the current quantum mechanical model¬ 
ing method. These discrepancies likely stem from many- 
body effects, such as exchange and correlation interac¬ 
tions (probably not correctly accounted for) and particu¬ 
larly excitonic effects (typically not regarded at all), when 
calculating the optical spectra of BiFeOs. 

The provided rigorous anisotropic optical modeling and 
the full dielectric function tensor for both BiFeOs phases 
will be highly useful to improve ab initio calculations to 
better understand multiferroic materials. 
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